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Abstract 

We present a general procedure for determining quasi-exact solvability of the Dirac 
^1 and the Pauli equation with an underlying sl{2) symmetry. This procedure makes full 

, use of the close connection between quasi-exactly solvable systems and supersymmetry. 

_ The Dirac-Pauli equation with spherical electric field is taken as an example to illustrate 

■ the procedure 
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1. In this talk we present a general procedure for determining quasi-exact solvability 

of the Dirac and the Pauli equation with an underlying sl{2) symmetry. This procedure 

I makes full use of the close connection between quasi-exactly solvable (QES) systems and 
H ■ 

5^ ! supersymmetry (SUSY), or equivalently, the factorizability of the equation. Based on this 
procedure, we have demonstrated that the Pauli and the Dirac equation coupled minimally 
with a vector potential [1, 2], neutral Dirac particles in external electric fields (which are 
equivalent to generalized Dirac oscillators) [3, 4], and Dirac equation with a Lorentz scalar 
potential [5] are physical examples of QES systems. 

Here we only give the main ideas of the procedures, and refer the readers to [1, 3, 5] for 
details. 



* Based on talks presented at the 11th International Conference on Symmetry Methods in Physics (Jun 
21-24, 2004, Prague) and at the XXIII International Conference on Differential Geometric Methods in 
Theoretical Physics (Aug 20-26, 2005, Nankai, Tianjin, China). 



2. For all the cases cited above, one can reduce the corresponding multi-component equa- 
tions to a set of one-variable equations possessing one-dimensional SUSY after separating 
the variables in a suitable coordinate system. Typically the set of equations takes the form 

+ = S-f., (2) 

where r is the basic variable, e.g. the radial coordinate, and f± are, say, the two components 
of the radial part of the Dirac wave function. The superpotential W is related to the external 
field configuration, and £^ involve the energy and mass of the particle. We can rewrite this 
set of equations as 

A-A+f_ = e/_, (3) 
A+A-U = eU, (4) 

with 

A^ = ±^ + W, e = S+£-. (5) 
ar 

Explicitly, the above equations read 

Here and below the prime means differentiation with respect to the basic variable. Eq.(6) 
clearly exhibits the SUSY structure of the system. The operators acting on f± in Eq.(6) are 
said to be factorizable, i.e. as products of A^ and y4+. The ground state, with e = 0, is 
given by one of the following two sets of equations: 

A+f^\r) = , /f (r)=0; (7) 
A-/f (r) = , f^^\r)^0, (8) 

depending on which solution is normalizable. 

One can determine the forms of the external field that admit exact solutions of the 
problem by comparing the forms of the superpotential W with those listed in Table (4.1) of 
[6]. 



Similarly, from Turbiner's classification of the sl{2) QES systems [7], one can determine 
the forms of W, and hence the forms of external fields admitting QES solutions based on 
sl{2) algebra. The main ideas of the procedures are outlined below. 



3. We shall concentrate only on solution of the upper component /_, which is assumed 

to have a normalizable zero energy state. 

Eq.(6) shows that /_ satisfies the Schrodinger equation = e/_, with 



with 



+ V{t) , (9) 



V{r) = W{rf - W'{r) . (10) 

We shall look for V{r) such that the system is QES. According to the theory of QES models, 
one first makes an "imaginary gauge transformation" on the function /_ 

/_(r) = 0(r)e-^W , (11) 

where g{r) is called the gauge function. The function 0(r) satisfies 



dr"^ dr 

For physical systems which we are interested in, the phase factor exp{—g{r)) is responsible for 
the asymptotic behaviors of the wave function so as to ensure normalizability. The function 
(j){r) satisfies a Schrodinger equation with a gauge transformed Hamiltonian 

Hg^-^ + ^Wo{r)-^ + [V{r) + VF^ - W^] , (13) 

where WQ{r) = g'{r). Now if V{r) is such that the quantal system is QES, that means 
the gauge transformed Hamiltonian Hq can be written as a quadratic combination of the 
generators J" of some Lie algebra with a finite dimensional representation. Within this finite 
dimensional Hilbert space the Hamiltonian Hq can be diagonalized, and therefore a finite 



number of eigenstates are solvable. For one-dimensional QES systems the most general Lie 
algebra is sl{2) . Hence if Eq.(13) is QES then it can be expressed as 

Hg^YI CabJ'^J^ + E C»>^" + constant , (14) 

where Cab, Ca are constant coefficients, and the J" are the generators of the Lie algebra 
sl{2) given by 

az 

f = iV = 0,l,2... 

Here the variables r and z are related by 2; = h{r), where h{-) is some (explicit or implicit) 
function . The value j — N/2 is called the weight of the differential representation of sl{2) 
algebra, and N is the degree of the eigenfunctions 0, which are polynomials in a {N + 1)- 
dimensional Hilbert space with the basis (1, 2;, 2;^, . . . , z^): 

(f) = {z - zi){z - Z2) ■ ■ -{z - zn) . (16) 

The requirement in Eq.(14) fixes V{r) and WQ{r)^ and Hq will have an algebraic sector 
with + 1 eigenvalues and eigenfunctions. For definiteness, we shall denote the potential V 
admitting + 1 QES states by V/v. From Eqs.(ll) and (16), the function /_ in this sector 
has the general form 

f-^{z- zi){z -Z2)---{z- zn) exp Wo{r)dr^ , (17) 

where Zi {i — 1,2, ... ,N) are N parameters that can be determined by plugging Eq.(16) 
into Eq.(12). The algebraic equations so obtained are called the Bethe ansatz equations 
corresponding to the QES problem [8, 1, 3] . Now one can rewrite Eq.(17) as 



/_ = exp(^- j Wj,{r,{z,})drj , (18) 

with 

iy^(r,{.,}) = iyo(r)-E^^- (19) 



There are possible functions Wjv(r, {zi}) for the sets of eigenfunctions 0. Inserting 
Eq.(18) into H_f_ = e/_, one sees that Wn satisfies the Ricatti equation 



where is the energy parameter corresponding to the eigenfunction /_ given in Eq.(17) for 
a particular set of N parameters {zi}. 

Prom Eqs.(9), (10) and (20) it is clear how one should proceed to determine the external 
fields so that the Dirac equation becomes QES based on sl{2): one needs only to determine 
the superpotentials W{r) according to Eq.(20) from the QES potentials V{r) classified in 
[7]. This is easily done by observing that the superpotential Wq corresponding to A?" = is 
related to the gauge function g{r) associated with a particular class of QES potential V{r) by 
g'{r) — Wo{r). This superpotential gives the field configuration that allows the weight zero 
{j — N — 0) state, i.e. the ground state, to be known in that class. The more interesting task 
is to obtain higher weight states (i.e. j > 0), which will include excited states. For weight 
j {N = 2j) states, this is achieved by forming the superpotential WN{r, {zi}) according to 
Eq.(19). Of the possible sets of solutions of the Bethe ansatz equations, the set of roots 
{zi, Z2, ■ ■ ■ , zn} to be used in Eq.(19) is chosen to be the set for which the energy parameter 
of the corresponding state is the lowest. 

4. Let us illustrate the above procedure by an example. We consider the motion of a 
neutral fermion of spin- 1/2 with mass m coupled non- minimally with an external electromag- 
netic field with an anomalous magnetic moment /i. The relevant equation describing such 
particle is the Dirac-Pauli equation [9]. This equation is useful in describing the celebrated 
Aharonov-Casher effect [10], and is also of some interest in quantum chromodynamics in 
connection with the problem of quark confinement [11]. 

We shall consider the situation in which only electric field E is present. In this case, the 
Dirac-Pauli equation Hip — Sip is described by the Hamiltonian 



W^-W^ = VN-eN , 



(20) 



H — ex. ■ p + ifx'-f ■ E + l3m 



(21) 
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with p = — iV and /3 = 7°. We choose the Dirac matrices in the standard representation 
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(22) 



where cr are the Pauh matrices. We also define -0 — {X: where t denotes transpose, and 
both X s-nd (p are two-component spinors. Then the Dirac-Pauh equation becomes 



<T • (p - ifJ,E)x ^ {S + m)ip , 
<r • (p + ii^E)(p = {S - m)x ■ 



(23) 



We now consider central electric field E = Eri. In this case, one can choose a complete 
set of observables to be (if, J^, J^, = 3/4, X). Here J is the total angular momentum 
J = L + S, where L is the orbital angular momentum, and S = |S is the spin operator. 
The operator K is defined a& K = 7°(5] • L + 1), which commutes with both H and J. 
Explicitly, we have 



K — diag {k^ —k^ , 
k = cr ■ L + 1 . 



(24) 



The common eigenstates can be written as 



here yj^n^O, 0) are the spin harmonics satisfying 



J 



and 



'^jyjmj 1 



k 



. _ 1 3 



(25) 



(26) 
(27) 
(28) 



(29) 
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where r is the unit radial vector. Eq.(23) then reduces to 

+ ^ + j /_ = {S + m)U, (30) 

(-|: + ^ + /^^.)/+ = {S-m)f.. (31) 

This shows that /_ and /+ forms a one-dimensional SUSY pairs with the supcrpotential W 
given by 

W^- + nEr , (32) 
r 

and the energy parameter e — £^ — m^. 

We can now classify the forms of the electric field Er (r) which allow exact and quasi-exact 
solutions. To be specific, we consider the situation where A; < and / dr/iEr > 0, so that 
/l^-* is normahzable, and /]°^ = 0. The other situation can be discussed similarly. In this 
case, Eq.(32) becomes 

W^-^-^ + I^Er. (33) 
r 

We determine the forms of Er that give exact/quasi-exact energy £ and the corresponding 
function /_. The corresponding function /+ is obtained using Eq.(30). 

5. Comparing the forms of the supcrpotential W in Eq.(33) with Table (4.1) in [6], one 
concludes that there are three forms of Er giving exact solutions of the problem : 

i) oscillator-like : iJiEr{r) oc r ; 

ii) Coulomb potential-like : iiEr{r) oc constant ; 

iii) zero field- like : iiEr{r) oc 1/r . 

Case (i) and (ii) had been considered in [12] and [13], and case (iii) in [12]. 
We mention here that the case with oscillator-like field, i.e. case (i), is none other than 
the spherical Dirac oscillator [11]. 

6. The form of the supcrpotential W in Eq.(33) fits into three classes, namely. Classes 
VII, VIII and IX of s/(2)-based QES systems in [7]. Below we shall illustrate our construction 
of QES electric fields in Class VII QES systems. 



The general potential in Class VII has the form 



yjv(r) = aV^ + 2abr^ + 6^ - a (4A^ + 27 + 3) + 7 (7 - 1) r"^ - 6 (27 + 1) , (34) 



where a, b and 7 are constants. The gauge function is 



q(r) — -r + -r — 7mr 



(35) 



We must have 0,7 > to ensure normalizability of the wave function. Eqs.(35) and (33), 
together with the relation Wo{r) = g'{r), give us the electric field E^: 



l_iE^°\r) = ar^ + br . 



(36) 



The Dirac-Pauli equation with this field configuration admits a QES ground state with 
energy S'^ — (e — 0) and ground state function /_ oc exp{—go{r)). Also, here we have 
1 — 1^1- 

To determine electric field configurations admitting QES potentials with higher 
weight, we need to obtain the Bethe ansatz equations for 4>. Letting z = h{r) = r^, Eq.(12) 
becomes 



-Az-^ + + 46^ - 2 (27 + 1)) - {AaNz + e) 



0(z) = . 



(37) 



For N — 0, the value of the e is e = 0. For higher N > and 0(r) = I\iLi{z — Zi) , the electric 
field El.^\r) is obtained from Eq.(19): 



/.Ef)(r)=/.EW(r)-X: ^'^'^ 



r[ h{r) - Zi ' 

For the present case, the roots Zj's are found from the Bethe ansatz equations 



(38) 



2az1 + 2bzi - (27 + 1) - 5^ 



0, i = l,...,7V, 



(39) 



and e in terms of the roots z/s is 



N 



e = 2(27+l)^-. 



i=i ■^i 



(40) 
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For = 1 the roots zi are 



and the values of e are 



e^^2(b± yjb^ + 2a(27 + 1)^ . (42) 

For a > 0, the root = < gives the ground state. With this root, one gets the 

superpotential 

WAr) = ar^ + br- ^ _, - - ■ (43) 
^ + 1% I 



From Eq.(38), the corresponding electric field is 



2r 



f^El^Ur) = ar^ + br - ^ , _, . (44) 

H + \zi I 

The QES potential appropriate for the problem is 

V{x) = Wl - W[ , 

= V^-e. (45) 

The one-dimensional SUSY sets the energy parameter of ground state at e = 0. Hence, 
the ground state and the excited state have energy parameter e = and e — — e' — 
4^62 + 2a(27 + l), and wave function 

/_ oc e-^°^'^ (r' - ^r) (46) 

and 

/_ (X e-^o^'') p - 4) , (47) 

respectively. 

QES potentials and electric fields for higher degree can be constructed in the same 
manner. 
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A more extensive discussion of the (quasi)-exact solvability of the Dirac equation in 
different background potentials can be found in [4]. 

This work was supported in part by the National Science Council of the Republic of China 
through Grant No. NSC 94-2112-M-032-007 and NSC 96-2112-M-032-007-MY3. 
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